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Abstract 

We re-examine the problem of strong coupling in a regularized 
version of DGP (or "brane-induced" ) gravity. We find that the 
regularization of ref. |hep-t h/0304148 differs from DGP in that it 
does not exhibit strong coupling or ghosts up to cubic order in the 
interactions. We suggest that the nonlocal nature of the theory, 
when written in terms of the 4-D metric, is a plausible reason for 
this phenomenon. Finally, we briefly discuss the possible behavior 
of the model at higher-order in perturbation theory. 
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1 Introduction 



Modifications of gravity in the UV or IR provide an interesting arena for ad- 
dressing phenomenological issues ranging from the hierarchy problem to the 
cosmological constant problem. In recent years, several ideas were proposed 
for a consistent modification of gravity in the infrared (pP-jS]). In particular, 
the DGP |2] model is a five-dimensional brane-world realization of a theory 
which gives modifications of four- dimensional gravity in the far IR. 

Higher codimensional generalizations of DGP were considered in |l]-[Zj; 
they differ from the co dimension-one model in several important aspects. 
For instance, they need UV regularization in order to properly define the 
propagation of gravity, since their Green's functions are UV divergent. Thus, 
it seems useful to study generalizations of these models which still share the 
same modifications of gravity in the IR, but are regularized in the UV. In jS] , 
we proposed a regularization procedure for the higher-codimension theories, 
which produces a ghost- and tachyon-free spectrum at linear order. The 
basic ingredient is to use purely higher dimensional gravity, with a region 
around the brane where the gravitational coupling is smaller (this may be 
dynamically realized by coupling gravity to a soliton; for discussions see 
refs.[2j, 8j- 10J). One can consider both smooth and sharp realizations of 
this mechanism; in this paper we shall confine ourselves to the sharp version, 
for which the action takes the form: 



with e a large dimensionless quantity and M* the fundamental Planck scale. 
One can add to this a brane term at the origin of transverse space, on which 
the Standard Model (SM) particles are confined. When one looks at the 
gravitational interaction between SM sources confined to the brane, using 
the Kaluza-Klein decomposition of the graviton, one sees that this action 
gives rise to four-dimensional behavior, within a critical distance scale and 
higher dimensional behavior above that scale [H]. For N=l this behavior 
coincides with DGP. 

In this paper, we further investigate this model in a somewhat different 
context. As was discussed in jTTj (see also [12]), the DGP model suffers 
from a problem also encountered in massive gravity, namely that the UV 
cutoff beyond which the effective field theory description breaks down is much 
lower than one would naively expect. In way reminiscent to how longitudinal 
polarizations of a spin one field with mass become strongly coupled at 
Ai = AirrriA/g, the longitudinal polarizations of a graviton with given mass 
m g become strongly coupled at a scale A 5 = (m^Mp) 1 / 5 . 
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A very useful method to find this scale in massive gravity using the 
Goldstone formalism was introduced in ^3]. In analogy with gauge theory, 
the authors identify the Goldstone bosons of broken coordinate invariance, 
which encode in a very simple way the scale at which the effective field the- 
ory description breaks down. These Goldstone modes are characterized in 
the Lagrangian formulation by fields with anomalously small kinetic terms 
(which, unlike the spin one case, vanish in the decoupling limit) and large 
self-interactions. For massive gravity, the Goldstone modes coincide with 
the longitudinal modes of the graviton at high energies. This is nothing but 
the gravitational version of the Goldstone equivalence theorem. In DGP, the 
Goldstone mode was identified to be the brane-bending mode which becomes 
strongly interacting at a scale A D qp = (Mp/L 2 ) 1 / 3 [TT], where L is the DGP 
scale; L = Mp/M^ . For typical DGP parameters (L of order of the Hubble 
scale), A^gp is around 1000 km, which is indeed unacceptable from the ex- 
perimental point of view. At this scale, from the effective field theory point 
of view, one loses control of the theory because one has to include all oper- 
ators consistent with the a priori symmetries of the model. This feature is 
therefore inherently a problem of quantum stability. 

At the classical level on the other hand, this problem translates into the 
appearance of the "Vainshtein scale" [Hj. Vainshtein noted that in mas- 
sive gravity, the linear approximation near a heavy source breaks down at a 
distance 

(i.2) 

For a black hole of mass M, for instance, this means that we cannot trust 
perturbation theory in massive gravity within d v , a distance much larger than 
the black hole's Schwarzschild radius (here we implicitly assumed the exper- 
imental constraint that the graviton mass is extremely small). This means 
physically that the nonlinear interactions associated with (massive) gravity 
are so large that, unlike the massless case, they can no longer be neglected 
at the scale d v . In the DGP model, as it was noted in ^Hj, the gravitational 
fluctuation around a heavy source contains a large term. Although this term 
is pure gauge at the linear level, it does contribute to the nonlinear correc- 
tions, which subsequently grow bigger than the linear fluctuations within a 
distance 

d DGP ~ (ML 2 /M|) 1/3 . (1.3) 

In this paper, we investigate how the above feature is changed in the regular- 
ized model. First, we review the main features of the model in the formalism 
of [TT], taking codimension N=l for simplicity. Then, we investigate how 
the model reacts to a heavy source located at the brane. We will see that, 
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in the regularized version, there is no blow-up mode in the linearized solu- 
tion, nor in the lowest-order (cubic) nonlinear interaction term, and that we 
can thus trust perturbation theory up until phenomenologically acceptable 
scales. The solutions do have the famous vDVZ discontinuity ^B], which 
means that the theory will give different predictions than standard Einstein 
gravity in the bending of light around the sun for instance, unless one intro- 
duces new dynamics to cure the discontinuity (using a repulsive vector force, 
for instance). 

Finally we go to the effective field theory description of all modes and look 
for modes which get strongly coupled at energies below the cutoff. We will 
see that the induced action on the brane is highly non-local in terms of the 
metric due to the non-trivial boundary conditions implied by the action (jl.lj) . 
Unlike the DGP case, all interacting modes receive healthy kinetic terms, 
signaling that the strong coupling problem effectively does not occur in our 
framework. So, our regularization has traded the strong coupling problem 
for the "almost locality" (in terms of the gravitational degrees of freedom) 
of DGP This is closely related to the behavior of "soft gravity" studied in 
ref We will comment on this and other features in the conclusion. 

2 A Scalar Toy Model 

Consider a scalar field described by the following action: 

S = -e f d b x [d M <pd M ip] - [ d 5 x [d M ipd M ip]. (2.4) 

J z<A Jz>A 

We are interested in the dynamics on the hypersurface z — 0. For that 
purpose we can calculate the effective four dimensional action for the scalar 
field ipo = <f(0) by integrating out the bulk [TT]. This procedure amounts to 
solving the bulk equations for ip with the constraint ip(0) = <po, and using 
this solution in the induced Lagrangian. Inside the "dielectric" (i.e. the 
region z < A), the solution to the bulk equations of motion is (we work in 
the four-dimensional Fourier transform for notational simplicity): 

(Pp(z) = cosh(pz)(fl + smh(pz)(p 2 p , (2.5) 

Outside the dielectric the solution which vanishes at infinity is: 

tf(z) (2.6) 

^DGP is local in the sense that at least the highest dimensional kinetic term operator 
is local in the gravitational fields. 



4 



We know that (p p is the boundary (four dimensional) field, which we keep 
fixed in our formulation. Solving the boundary condition at z = A, we get: 



2 1 + etanh(pA) ^ 
^ = " e + tanh(pA) ^ 

Now the induced action at z = is the boundary term obtained by partial 
integration of equation (j'2.4j) : 

Sfewfa) = t l d^X [lf p (z)d z lf p {z)]\ . (2.8) 



Since d z (p p (z) = pip p and by going back to position space we obtain the 
following induced action: 

"1 + etanh(v/^n7A) 
_ e + tanh( v /^D7A) 



Sbrane{<P) = ~£ J VW- Q 4 



V9 . (2.9) 



This action looks highly non-local, but we can use it as an effective field 
theory with cutoff Kuv — A -1 . Then we see that for pA <C 1 the action is: 



(2.10) 



So what is the net effect of the five dimensional dielectric on the four dimen- 
sional induced action? There appear two regimes on the brane, defined by the 
"crossover scale" L = eA. In the regime pL 3> 1, there is four dimensional 
propagation with propagators scaling as 1/p 2 , whereas in the complementary 
regime pL 1, the propagation is five dimensional. So by this mechanism 
we get the same propagation as in the DGP model. The crucial difference 
here is that we do not have to include a separate four dimensional kinetic 
term on the brane to get four dimensional behavior within a certain range 
at the brane. It is induced by the bulk Lagrangian itself. 



3 The Vainshtein Scale 



We look at the gravitational action on a half space, in which there are two 
bulk regions to be distinguished: inside and outside the "dielectric". The 
gauge-fixed bulk action reads: 



Ssuik/Ml — eS z< & + S z> a, 



where SWa is given by: 



S x 



R(j) - K^KV + K 2 — -F* 



(3.11) 



(3.12) 
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Here K^ v is the extrinsic curvature and we used the notation = and 
N = # 5 ~ 5 1/2 . The mass scale M* is the bulk Planck scale (~ 10- 3 eV-100MeV), 
and e = Mp/M^. For A we choose the natural value M~ x . Finally, we 
choose the de Donder gauge fixing in the bulk F M = djsi{Jgg MN ). At the 
linearized level, expanding g MN = r] MN + Hmn, this action gives the usual 
bulk equations for the fluctuations: 



a 5 (h 



MN 



~VMNh) 



0. 



(3.13) 



This equation holds inside and outside the dielectric, and the two solutions 
of this equation should be matched over the surface using the appropriate 
boundary conditions. Varying the action in the two regions gives separate 
contributions to the interface z = A, which give us the following conditions 
[we use the notation [J\ = eJ(A_) — J(A + )]: 



[K^ u - % V K - -V^uF 5 \ = 0, 



(3.14) 



|F M J=0, (3.15) 

[F 5 \=0. (3.16) 

These equations are obtained by variation of h^, and iV respectively. F M 
appears here in its familiar linearized form; F N = du{h MN — \r] MN h^). Using 
for instance the trace of equations ()3.14|) and JHHEJ), we obtain (x = d z x): 



VJ = 0. 

On the other hand, equation ()3.14|) gives: 

V - 2d {ll N v) \ = 0. 



(3.17) 



(3.18) 



If we assume a source which is localized at the brane, the boundary conditions 
at the brane are given by the following set of equations: 



hpv - 2d {lM N u) - T)^ ( -/i 4 - d.N + -h 5 5 



1 T, 



2 eAf 3 ' 



d K V + N„ - -8^ {h 4 + h 55 ) 



0. 



(3.19) 

(3.20) 
(3.21) 



These are the initial conditions on the set of gravitational fields. Now we no- 
tice that the bulk linearized equations of motion for and h 55 , C^iV^ = 
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and □ 5 (/i 55 — h^) = respectively, allow us to make the following consis- 
tent gauge choice; = and h 55 = /i 4 (this is in specific consistent with 
the boundary conditions given by equations Ij3.19j) - (|3.21j) ). This means we 
can essentially treat this as a scalar system at the linearized level, with in 
particular the fluctuation obeying: 



LVJ = o. 



(3.22) 



Assume now that we couple gravity to a source on the brane. What would 
the solution for look like then? The boundary condition at the brane is: 



(3.23) 



This leads to a solution of the form (in four-dimensional momentum space): 



h^u{p,z) = --D(p,z) 



(3.24) 



Here D(p, z) is the bulk to boundary scalar propagator of the theory, which 
is calculated in the appendix. The point is made by contrasting the solution 
obtained in the DGP model, using the de Donder gauge in the bulk: 



hau(p, z) = --D'(p,z) 



t ^{p) ~ ^r)^T(p) - ^p^p u T(p) 



(3.25) 



where D'(p, z) 



'JT^p 2 +p/L 



e pz . The difference is clear: in the DGP case we 



get an extra term which blows up rapidly with increasing p, while this term 
is absent in the regularized model. For DGP, this mode does not contribute 
at linear order, since it is pure gauge, but it does contribute at quadratic 
order (cubic terms in the Einstein-Hilbert action). This means that in DGP, 
as in massive gravity, the linearized approximation breaks down at a distance 
scale much higher than the massless case. In DGP and massive gravity, the 
contribution from the gravitational radiation to the energy density is much 
bigger than in the massless case, so that the heavy back-reaction in these 
models is not negligible. We discuss this in more detail now. 

A source confined to the brane, Tmn(x, z) = 6(z)6 l ^[8^ r T fll ,(x) gives rise at 
linear level to h^(T), given by equation ()3.24j) . This fluctuation contributes 
now by itself to the energy-momentum tensor T^ N (x, z) = /i < - 1 ^(T)G/i^ 1 ^ ) (T), 
with G some operator with complex tensor structure but quadratic in the 
derivatives. This gravitational radiation is now everywhere in space, not just 
confined to the brane. Perturbatively we can now solve for h^ 2 \ which will 
be quadratic in the source. There are obviously two possibilities. The first 
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one, like in standard Einstein gravity, is that is very small in comparison 
to so that perturbation theory can be trusted, and one can keep the 
linearized approximation up until, say, the Schwarzschild radius of a black 
hole. The second possibility is that grows bigger than at some large 
distance scale. This can be due either to some large pieces in which are 
pure gauge only at linear order [and can be present in physical contributions 
to h^], or because the propagator of the system contains large pieces, which, 
applied to T^ N generate fast growing terms. 

Instead of dwelling on this point, we may as well compute the second order 
fluctuations and see if either one of the possibilities is present. To solve for 
this system in the de Donder gauge we need the scalar propagator, which is 
calculated in the appendix. The equations of motion take the following form: 



a 5 (h 



(2) 
MN 



-riMNhf ] ) = 



^iT^ N (x, z), 



(3.26) 



where the coupling ki is different on both sides of the interface. Obviously 
the extra constraints = and h^j = hf^ are not automatically satisfied, 
but we can use residual gauge invariance to put: 



V i 

(2) I 



o, 

(2) I 
55 |A- 



(3.27) 

h^'\A = h^\ A . (3.28) 

Indeed, any coordinate transformation £ M (x,z) obeying 5 £ M (x,z) = will 
leave the equations of motion invariant (but not the boundary conditions), 
and they can be used to enforce the above identities. Since 5Hmn = 2<9(m£jv) 
is a gauge transformation, it will not contribute to any amplitude between 
sources since we assume them conserved. We will discuss these transforma- 
tions in more detail below. In any case, we see that the boundary conditions 
for h$ and hf§ are now the same as before, so that the scalar propagator 
still determines the solution uniquely: 



h 



(2) 
MN 



^ J d b x'D(x-x';z,z') 



T$ N (x',z')-lr]MNTM(x',z') 



(3.29) 



In other words, if neither the scalar propagator D(x — x'; z, z') (see appendix) 
nor the source grow big anywhere, the second order solution does not grow 
big anywhere. So we can be sure that the linearized approximation does not 
break down at any unacceptably low scale. 

In the preceding derivation, we have assumed that the gauge transforma- 
tions we use are globally defined, i.e. that they are proper gauge transforma- 
tions. One might suspect that in fact the interface will bring complication, 



8 



and one has to check that the gauge degrees of freedom truly decouple globally 
(from T^). This amounts to checking if the following integral vanishes: 

S = [ d 5 x & M t^Tg> N + [ d 5 x d^^TSi- (3.30) 

J z<A Jz>A 

Since gravitational radiation is locally conserved: 

d M T$ N = 0, (3.31) 
the only contributions left are exactly the ones at the interface: 

S, = f d A x eT$\ A - - J d 4 x t*T$\ A+ , (3.32) 



S* = 



J d 4 x £ 5 T 5 (2) | A _ - J d'x e 5 T 5 ( 5 2) | A+ . (3.33) 

In the first integral, we can calculate the directly. Since = 

everywhere, the expression for is particularly simple. The reader can 
check that, indeed: 

T^-R^l (3.34) 

where is the part of the Ricci tensor quadratic in h$. The constant of 
proportionality is the inverse coupling constant: eAlf inside the dielectric, 
and Ml outside of it. Explicitly: 

<s = ~\h Xv d^h Xv + l -h x »d b d x h, v + ^d x h 55 d 5 h* 

-\d,h MN 8,h MN . (3.35) 
Now we notice that all terms appearing are linear in d 5 h,MN, which means: 
T$\ A _ = eM^gU- = MlR%\ A+ = T$\ A+ . (3.36) 

This was expected if our system was to be stable, of course: any inequality 
here would imply nonzero pressure on either side of the dielectric. Since we 
use continuous gauge transformations we can conclude that = 0. 

As for S5, we only have to notice that £5 vanishes at the boundary. Indeed, 
since we demanded that /155 is continuous over the boundary and [^55 J = 0, 
we can conclude that the gauge parameter £ 5 (which transforms d~h$ 5 = 2<9 5 £ 5 ) 
obeys [£5 J = 0, implying £ 5 | A = 0. 

The vanishing of S5 and also ensures that, at cubic order in the inter- 
actions, all the ghost degrees of freedom decouple. It would be interesting to 
see if this property persists to all order. 
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4 Mode Analysis 



We now turn to the mode analysis of the complete four-dimensional induced 
Lagrangian. Before we start, it is a good idea to review shortly what happens 
in the case of massive gravity and DGP. 

In massive gravity, the only local, ghost-free, free Lagrangian is given by 
adding the Pauli-Fierz term to the linearized Einstein-Hilbert action: 

C PF = M 2 ^R L (g)\ g=v+h + m 2 g M 2 P (hl v - h 2 ). (4.37) 

The Pauli-Fierz term breaks coordinate invariance explicitly, and this La- 
grangian propagates five degrees of freedom. One can restore this invariance 
by introducing a Stuckelberg field, by the substitution h^ u — > h^ v + 2d( p A u ), 
and demanding that under coordinate transformations SA^ = — £ M ( [TT| . [T5] ) . 
We decompose A^ in a longitudinal and transverse part A^ = A^ + d^a. 
A^ is essentially the Goldstone field of broken coordinate invariance, whose 
proper nonlinear definition was given in [13] ■ This substitution leaves the 
Einstein-Hilbert action invariant, as it is invariant under coordinate trans- 
formations, but in the mass term A^ appears in the following form: 

C mass = m 2 g M 2 [F 2 u + 2A^h? v - 2A^h v l/ ], (4.38) 

where is the field strength of A^. One now immediately sees that the 
transverse part of the field has the correct kinetic term, but that the longi- 
tudinal scalar gets a kinetic term only by mixing with h^. The resulting 
kinetic term, after diagonalization, reads: 

L a = m A g M 2 P aU A a. (4.39) 

Now the logic is as follows. The field a is expected to obtain its kinetic 
term only by mixing because otherwise it would obtain a higher derivative 
kinetic term, implying violation of unitarity. The kinetic term it obtains is 
anomalously small as a consequence of the mixing, which means that the 
higher dimensional operators containing a, which are generated in the inter- 
acting theory, could actually become important at a very low scale. To see 
when these interactions become important we have to use the canonically 
normalized field a c = a/(m 2 Mp). At the quantum level, the only symmetry 
protecting operators being created is the shift symmetry mentioned above. 
Thus we expect for instance the appearance of interactions of the form: 




(4.40) 
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This term can also be seen to appear explicitly in the self-interaction terms 
of the Goldstone boson ^3]. This interaction is in fact the strongest in- 
teraction that will appear in the Lagrangian, and it becomes strong at 
A 5 = (mjMp) 1 / 5 . At this scale we should include all operators consistent 
with the symmetries, and we need a UV completion to regain predictability 
above that scale. 

For DGP, the situation looks vastly different at first. The theory obeys 
manifestly general covariance, and hence it does not share the arbitrariness 
of massive gravity. Recall that the Pauli-Fierz term is only uniquely de- 
fined at quadratic order in the fluctuations, and its nonlinear extensions are 
quite arbitrary. Specifically, the generation of quantum interactions beyond 
quadratic order is not protected by general covariance. The induced gauge- 
fixed action at the brane takes the form [TT] : 



£>dgp — Mp 



L 2^( n 4 - mA) V - \h 4 ( D 4 - mA)h 4 



in Tft 
-mJV"(A + m)N^ + -h 4 Ah 55 - mN»d^h 55 - -/i 55 A/i 55 J , (4.41) 

where we used the notation A = y/— D 4 , and where m — 1/L — M^/M P . 
By analyzing the scalar modes in this Lagrangian , one sees that there is 
indeed a scalar mode it, defined by = d^n, h 55 = — 2 An 1 1J , which gets 
a kinetic term only by mixing with the trace part of h^ u . The mixing term 
in the action f)4.41j) is proportional to: 

M|m7rn 4 /i 4 . (4.42) 

This means that after diagonalization this mode will receive a kinetic term: 

~ Mj,m 2 7rD 4 7r. (4.43) 

Again we see that the kinetic term is suppressed. Now, the bulk interactions 
induce the following cubic interaction at the brane: 

C int = Mpmd^nd^U^. (4.44) 

Using the canonically normalized field, ir c = 7i/(Mpm), one can easily com- 
pute the scale at which this interaction becomes important: A3 = (Mp/L 2 ) 1 ^ 3 . 
Unlike in massive gravity, this operator has canonical dimension seven, it is 
explicitly present in the classical Lagrangian and is not renormalized. The 
mode has, by its definition, the physical interpretation of being the "brane 
bending." 

We are now in the position to investigate the appearance of a strong 
coupling problem by following the procedure used in the previous examples. 
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We will calculate the induced Lagrangian of our system and investigate all 
its modes. We may expect some scalars to receive their kinetic terms by 
mixing, especially the brane bending mode. If this mixing is indeed small, 
we can conclude that our model shares the same problem as massive gravity 
and DGP. 

In order to calculate the induced Lagrangian, we need the on-shell bulk 
solutions of each mode (in the de Donder gauge), with all the fields con- 
strained by the boundary conditions given by equations 1)3. 14113. IB)) . Using 
the orthogonal decomposition of the four dimensional fluctuation into its 
transverse-traceless, longitudinal, and trace parts 

V = hT J + 2 ^C) + V»v<f>, (4-45) 
we get the following (mixed) boundary condition from equation (j3.14|) : 

L^J + W>- 2 ^)J =0, (4.46) 

where we defined for later convenience i]^ = N^ — (^. Taking the trace of this 
and comparing with the double (four dimensional) divergence of equation 
(j3.18|) we see that: 

L0J = 0. (4.47) 

Now it is easy to obtain the remaining boundary conditions on the fields. 
For the transverse traceless part we have: 

[h T J\ = 0. (4.48) 

For Nfj, and we essentially get conditions on their linear combination 77^: 

M = 0, (4.49) 

lVn ~ = • (4-50) 

Here, we kept only the leading term in an 6 > 1 expansion. So, we notice 
that the transverse and trace part of h^ u as well as h 55 obey homogeneous 
boundary conditions, while iV M = NT + 8^9 and obey mixed boundary 
conditions. Furthermore, we ask continuity of the fields h^, h 55 and iV M 
(this last constraint follows from imposing the gauge fixing). For a field 
obeying homogeneous boundary conditions (continuity of J and [^J = 0), 
the solution is the same as the scalar example discussed in sections 2 and 3. 

§This decomposition is unique because 4-D fields that are transverse and longitudinal 
at the same time have vanishing 4-D momentum. So, they do not extend to 5-D on-shell 
fields vanishing at z = 00. 
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Those fields thus include the transverse traceless part of h^ u , h^J, and two 
scalars, and h 55 . We proceed by solving for 77. It obeys equation ()4.49|) and 
has a continuous derivative over the interface. Setting r/^ = rj^ + S^ry, we see 
that only the longitudinal part obeys non-homogeneous boundary conditions, 
and we can solve it by: 

rj = rjj sinhQoz + b) z < A, (4-51) 
v = r] e- pz z > A. (4.52) 

We easily find that b = — 1/e — pA. By equation (|4.5()|) and the continuity 
of 4>, h 55 and 77, we see that these fields are related by: 

h 55 \ « /i 55 |a « 2(77 - 0)| A « 2(r) - 0)| o . (4.53) 

This constraint already signals a major difference between our regularization 
and DGP in the de Donder gauge. In DGP, the boundary theory has three 
independent scalar fields, while here there are only two [see eq. ( 14.41)1 and 
ref. [□]]. 

The induced brane action obtained by integrating out the bulk linearized 
action ()3.11|) reads: 



ind j 



d x 



- h^hfu, - -h 55 h 5 5 + ^h 4 h 4 - 2N fl N> M 



+l(h 55 h 4 + h 4 h 55 ) - 2N^2d K h K » - d»h 4 - d»h 55 ) ( 1.5 j) 







The transverse traceless part and the transverse vector in the metric decouple 

TT 



from the other modes. Concentrating on hT£ for instance, it will appear in 



the action only with the term: 

T md (h TT ) = f d'x [h TT ^(B 4 - A/L)h T J] . (4.55) 

We further concentrate on the scalar sector of the theory. Substituting all 
the projections of the field, one obtains after some algebraic manipulations: 

fM 3 f .... 
IWC, 0, ^55, 9) = f- J d*x [-2CQ 2 C + 20DC + 20DC + 4# 

-7^55^55 + h 55 D( + h 55 n ( + 2h 55 df) + 2h 55 (fi + 

2» + 4#n 2 C-4#D0-2#n/i 55 ]. (4.56) 

We now define: 

if, = rj = + n 4 C, (4.57) 
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where the last equation follows from the 5-D on-shell condition. We can 
eliminate /155 from the action using the constraint (|4.53J1 . so that the action 
gets the simple form: 

eM 3 f ... 
0) = f- \ dx [- 6 # + W + 6^0]- (4.58) 

We notice that the field ip does not appear to have a diagonal kinetic term. 
Since the fields ip and <ft have the same behavior [in momentum space ~ 
Acosh(— 1/e — pA +pz)], the last two terms are identical, and we are ready 
to write the final form of the action using the solution of both fields in the 
bulk; 

IW^, 4>) = \J d 4 x [3M|0(D 4 - mA)0 + 6eM^n 4 # (4.59) 

Indeed, the scalar 7] is the one we set out to find from the beginning. As in 
both DGP and massive gravity, it gets its kinetic term by mixing with the 
scalar part of the metric. But unlike in DGP or massive gravity, the mixing 
term is very large] It will not give rise to an anomalously small kinetic term, 
but a very big kinetic term for the field rj; 

Cr, ~ M^M^7]D A r]. (4.60) 

The canonically normalized field r] c = r]/(MpM*) will not cause trilinear 
interactions in the theory to blow up at unacceptably low energy scales. As 
a rough estimate, we can see that the most dangerous induced interaction 
will give a cubic coupling of the form: 

£cub ~ M p M 2 VcP 2 VcP 2 Vc (4.61) 
which gives a strong coupling scale: 

A 3 = (MpM, 2 ) 173 > M*. (4.62) 
So we can trust our framework up until the cutoff we imposed, M*. 



5 Conclusion 

The results obtained in this paper are puzzling. We could have thought that 
our regularization of the DGP model differs from the unregulated version 
only at distances D « A. Instead, we have seen that it differs even at 
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the vastly larger scale A^ GP . How is this possible, and, where are potential 
problems hidden? 

The origin of this difference can be seen in two ways. First, by comparing 
eq. |L21 with the DGP result eq. In eq. (jS21, unlike in the DGP 

case, the trace of the metric fluctuation, M, propagates and is a ghost. This 
is not a disaster because the difference between the DGP propagator and that 
of the regularized theory is a gauge transformation. The absence of a gauge 
mode, that becomes large at the scale d v , is also a feature of the "softly mas- 
sive" gravity studied in ref. So, as long as linear gauge transformations 
decouple, all ghosts decouple. In section 3, we explicitly checked that this 
indeed happens to cubic order in the interactions. Clearly, an all-order check 
of this property is crucial for the ultimate viability of our regularization. 

The other way of understanding the difference between DGP and the 
model presented here is to notice that the linearized action eqs. (|4.55|4.5(jjl 
is very nonlocal when expressed in terms of h^ v . The free DGP action is 
local in the UV, since its highest-dimension term is local. The action in 
eqs. (I4.55l4.56p . instead, is never local because of the implicit presence of 
spin-projection operators. By this we mean the following: a high-spin theory 
is always nonlocal, when written in terms of just its physical degrees of 
freedom. In massive gravity or DGP, one can add unphysical polarizations 
-a massless vector and a scalar for each helicity ±2 polarization- so that the 
free kinetic term is local. In our case, this is not possible, because of the 
uniqueness of the Pauli-Fierz action. 

When written in terms of the spin-projected components, the action still 
exhibits a potential problem due to the presence of ghosts. It is sufficient 
to notice that in the free action eq. fj4.58j1 . the field x = $ — CUC never 
appears. It has zero kinetic term, so, if it appears at all in the nonlinear 
terms, it would interact with infinite strength. Since x is precisely a 4-D 
gauge transformation, our calculation in section 3 proves again that this does 
not happen at cubic order in the interactions. This gauge transformation only 
acts on the 4-D coordinates, while leaving the interface at z = A fixed. So, 
it is reasonable to think that it is non-anomalous, and that it will decouple 
at all orders. Again, an explicit proof of this property is important for the 
consistency of our regularization. We hope to address this question in the 
near future. 
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Appendix A: The Scalar Propagator 



In this section, we calculate the scalar propagator for our system. We basi- 
cally have to solve the following equation: 

tf-%)4> = U(z-*>). (5.63) 

There will be two distinct cases: z' < A, hi = eMf , and z' > A, hi = M%. 
We treat them separately here. 



Case I: z' < A 



There are three different regions; region I, which is for values < z < z' 
(we demand that 0|o = 0), region II, which is for values z' < z < A, and 
then region III, which is the space outside the dielectric. The solutions to 
equation (|5.6H|) are: 



(f>i(p,z) = Acosh(pz), 
4>u{p, z) = B cosh(pz) + C smh(pz), 
<f> III (p,z) = De-*( g - A >. 

The continuity equations are: 

A cosh(pz') = B cosh(p/) + C sinh(pz'), 
D = £cosh(pA) + Csinh(pA), 

while the derivative boundary conditions read: 

pB smh(pz') + pC cosh(j»2;') — pA sinh(pz') 



(5.64) 



(5.65) 



(5.66) 



e [pB sinh(pA) + pC cosh(pA)] = —Dp. 



Adding the continuity and derivative conditions we obtain, in the regime 
where pA <C 1 and in the approximation e ^> 1: 



Acosh(pz') = B 



cosh(py') — 



C 

1 

e 



- + Ap)B, 



Ap ] sinh.(pz') 



(5.67) 



B [psinh(pz') — p 2 A cosh(p2')] — Apsiah.{jpz') 
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The last two equations give us A ps B and: 

1 1 



eAM 3 p 2 + jl> 

c = 1 



(5.68) 



eMl V 

This leads to the following propagator for z' < A: 

(f>i(p, z) = (pii(p, z) = - ^ 2 ^ 2 _^ p / L cosh(pz), (5.69) 

Case II: z' > A 

There are three different regions; region I, which is for values < z < A (we 
demand that <p\o = 0), region II, which is for values A < z < z', and region 
III, for which z' < z. Now the solutions are as follows: 

(f>i(p,z) = Acosh(pz), 
4>ii{Pi z) = B cos\i{pz) + C sinh(pz), (5.70) 
<j )lII (p,z)=De-^- z '\ 

The continuity and derivative equations read: 

Acosh(pA) = Scosh(pA) + Csinh(pA), (5.71) 
D = B cosh(pz') + C sinh(pz'); 



epAsinh(pA) = psmh{p/S.)B + pcosh(pA)C, (5.72) 
—Dp — pB sinh^z') — pC cosh(pz') = — 

Adding the derivative and continuity equations leads after some simplifica- 
tions to A ps and: 

C ps eApB, (5.73) 

5 = -- - e~ pz ' 

Ml p + p 2 L 
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We write down the explicit solution for the propagator at z' > A: 

1 1 / 

z ) = -j^^^z e ~ pz ' cosh (^)> 

lli 1 pL i 

(f>n(p,z) = - e~ pz ' cosh(pz) - e~ vz ' sinh(pz), (5.74) 

M% p + p z L M£ p + p z L 

(f>in(p,z) = - 775 —07 [cosh(pz') + sinh(pz')] e _p2 . 
M£ p + p z L 

We see that in the case pz' ^> 1: 

fe(^) = -^ pM) . (5.75) 



The bulk to boundary propagator is obtained by putting z = 0: 

1 1 



>i = 911 



M 2 P p 2 +p/L : 

i 1 1 r -W 

m M 2 p p z +p/L 



(5.76) 



The propagator is well-behaved everywhere and does not grow anomalously 
large in any region. 
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